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ABSTRACT 



Aims. We present a study of the envelope morphology of the carbon Mira R For with VLTI/MIDI. This object is one of the few 
asymptotic giant branch (AGB) stars that underwent a dust-obscuration event. The cause of such events is still a matter of discussion. 
Several symmetric and asymmetric scenarios have been suggested in the literature. 

Methods. Mid-infrared interferometric observations were obtained separated by two years. The observations probe different depths 
of the atmosphere and cover different pulsation phases. The visibilities and the differential phases were interpreted using GEM-FIND, 
a tool for fitting spectrally dispersed interferometric observations with the help of wavelength-dependent geometric models. 



Results. We report the detection of an asymmetric structure revealed through the MIDI differential phase. This asymmetry is observed 
at the same baseline and position angle two years later. The observations are best simulated with a model that includes a uniform-disc 
plus a Gaussian envelope plus a point-source. The geometric model can reproduce both the visibilities and the differential phase 
signatures. 

Conclusions. Our MIDI data favour explanations of the R For obscuration event that are based on an asymmetric geometry. We 
clearly detect a photocentre shift between the star and the strongly resolved dust component. This might be caused by a dust clump 
or a substellar companion. However, the available observations do not allow us to distinguish between the two options. The finding 
has strong implications for future studies of the geometry of the envelope of AGB stars: if this is a binary, are all AGB stars that show 
an obscuration event binaries as well? Or are we looking at asymmetric mass-loss processes (i.e. dusty clumps) in the inner part of a 
carbon-rich Mira? 

>: 
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■ 1 . Introduction 1. the ejection of a spherical shell by the star; 

2. enhanced dust-condensation in an already existing spherical 

] The C-star R For is classified as a Mira variable in the General shell' 

■ Catalogue of Variable Stars ( |Samus et al. |,|200§. Objects of this 3 an as ' ymme tric dust ejection in a preferential direction (i.e. a 
variability class show photometric variations with amplitudes ( jj sc or a more complex structure)- 
larger than 2.5 mag in the V band, and with long periods (> 100 4 an a ' symme tric dust ejection in a random direction, 
days). The light curves of Mira variables are usually described as 

$J ; regular with sinusoidal behaviour. Nevertheless, long-term ob- | Le Bertrd (fl988h and IWinters et al.l (fl994l) supported the first 
servations of theseobjects show erratic drops in their brightness two scenarios, which reproduced the photometrical data with 
( iBarthes et all [19961 and references therein). dynamic models in spherica l symmetry. Advan ced 2D and 
The var i ability of R For was discussed for the first time by 3D model atmospheres by IWoitke & Niccolinil (|2005) and 
|Feast et alj (|1984|), and later by |Le Bertre| (|1988|). Both works |Frevtag & Hofnerl d2008l) predicted the formation of concen- 
concluded that the light curve showed a deep minimum in 1983. trated dust clouds in the inner part of the atmospheres, and a 
This decrease was attributed to an increased dust obscuration, spherical distribution on a large scale. Such a morphology might 
the reason for which is (still) unknown. explain the dust obscuration events. However, up to now those 
Four scenarios were envisaged by iFeast et all d 19841) to ex- models were not compared with observations. On the other hand, 
plain the obscuration event, the first two scenarios were symmet- these obscuration events have simil arities with those ob served 
ric, and the latter two involved an asymmetric geometry of the for R Coronae Borealis stars (R CrB; IClavton et al.Ll2012l) . In the 
envelope. The scenarios involve: latter case, th e third and the fourth scenario would be more likely 
explanations (Bright et all 120 1 lUJeffers et al.Ll2012h.lLeao et al.l 
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* Based on observations made with ESO telescopes at La Silla (120071) used VLTI/MIDI to study the circumstellar environment 
Paranal Observatory under program IDs 080.D-0231 and 084.D0361. of the R CrB star RY Sgr. This star showed a drop in the J 
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and H light curve similar to the one observed for R For. The 
observations of RY Sgr showed that the star is surrounded by 
a single dusty cloud l ocated at ~100 stellar radii (~30 AU). 
Whiteloc k et alJ d 19971) suggested a binary-related effect as an 
alternative explanation. Exceptionally deep minima in the light 
curves of single O-rich Miras are not observed, but they are fairly 
common for symbiotic Miras. 

This letter presents high angular resolution observations in 
the mid-infrared that allow us to probe the dust-forming region 
of the nearby AGB star R For. Observations close in time, with 
different position angles and baseline lengths, the differential 
phase information, and the comparison with geometric mod- 
els can help to distinguish possible deviations from a spherical 
structure. This allows one to confirm or reject the different sce- 
narios suggested in the past. Observations and data reduction 
procedures are presented in Sect. [2] The study of the morphol- 
ogy by means of geometric models is described in Sect. [3] and 
is discussed in Sect. [4] 

2. Observations and data reduction 

To study the morphology of the close circumstellar environment 
of R For, high angular resolution measurements were acquired 
in the mid-infrared with the ESO-VLTI/MIDI facility (program 
IDs 080.D-0231 and 084.D0361, Ohnaka and Sacuto respec- 
tively). In Table Q] the 29 visibility points observed for the sci- 
ence target are presented. Twenty-four points out of the original 
sample turned out to be of good quality. In this work we con- 
centrate on the data marked in Table Q] with (b) and (d): these 
observations were obtained (within each group) at roughly the 
same pulsation phase (+ 0.03), but they sample different spa- 
tial scales of the object. Therefore, these two groups are ideal 
for studying the morphology of the star. The remaining obser- 
vations will be presented in a forthcoming publication (Paladini 
et al., in pre p.), where we compar e them with dynamic model 
atmospheres (Matt sson et alllToiOl) . 

The uv coverage dispersed in wavelength used for this work 
is plotted in Fig.Q] Every observation of R For is paired with the 
observation of at least one calibrator. The list of calibrators used 
for the data reduction is given in Table [2] The data were reduced 
with the standard data reduction pipelines MIDI Interactive 
Analysis ( MIA), and Expert Workstat ion (EWS) package ver- 
sion 1.7.1 (iKohlen. l2005t iJaffel |2004). A detailed description 
of the data quality check and error determination is given in 
lKlotzetalJ(l20l2al) . 

The derived visibilities (upper panels of Fig. |2]i exhibit the 
typical tilde-shape of A^-band interferometric o bservations of 
C-stars with a SiC featu re present at 11.3 fim dOhnaka et al.L 
l2007t ISacuto et all l20lTh . We recognise a drop in the 8 - 9/im 
region, where the main contributors to the opacity are C2H2 and 
HCN (Fig. I}. We note that the long baseline visibilities of 2009 
are systematically lower than the observations of 2007, indicat- 
ing a more extended dust shell in 2009. The differential phase 
is a crucial observable for this work because it carries informa- 
tion about the morphology of the source. This corresponds to the 
difference in phase between the different spectral channels. The 
EWS software performs a linear fit to the index of water vapour 
refraction to clean the phase s ignal from the atmospheric effects 
(IJaffel l2004fc frubbs et all l200l . In R For all differential phases 
show zero signal, except for those observed with the longest 
baselines (see lower panels of Fig. 0. The phase information 
in these observations clearly points to a non-central symmetric 
structure. Signatures in the differential phase have been reported 
in the literature for two silicate J-type C-stars: IRAS 18006-3213 



bv iDeroo et al ] (120071) . and BM Gem bv lOhnaka etaD (12008). 
These objects show silicate emission that probably originates 
from a circum-binary or a circum-companion dis c. Many C-stars 
are known to show photosph eric asymmetries dRagland et all 
12006b ICruzalebes et all l2012h . while very few extreme objects 
depart from spherical m orphology in the inner dust shell (e.g . 
IRC+10216 and CIT6: McCarthy et all 1 19801: IWeieelt et all 
1998; iMonnier et all 120001: iTuthill et all 120001: Ichandler et ail 
2007). To our knowledge this is the first time that MIDI has de- 
tected a signature of deviation from central symmetry for a car- 
bon Mira variable through the differential phase. 

3. Morphological interpretation 

The w avelength-depende nt geometrical model-fitting tool GEM- 
FIND ( Kl otz et alll2012bl) was used to model the interferometric 
observations obtained at comparable visual phases in 2009 ((d) 
in Table Q]). GEM-FIND is based on the Levenberg-Marquardt 
least-squares minimisation method and allows the fitting of dif- 
ferent centro-symmetric and asymmetric geometric models to 
the dispersed calibrated visibilities. The fitting is limited to the 
wavelength range < 12.5 fim to avoid any bias that may be 
caused by the noise that usually affects MIDI measurements at 
longer wavelengths. The best-fitting model visibilities were then 
used to derive synthetic differential phases, which were com- 
pared to the observed ones. 

Spherical and aspherical optically thick (i.e one-component) 
and optically thin (i.e., two-component) models were tested. The 
best-fitting model (Fig. [3) consists of (i) a uniform disc (UD) that 
represents the central star, and (ii) an extended centrosymmetric 
Gaussian distribution that represents the envelope. Furthermore, 
to account for the non-zero differential phases at long baselines, 
(iii) a third component corresponding to a Dirac delta function 
was added. This geometric model is defined by six parameters 
listed in Table [3] with the resulting best-fitting values (or range 
of values for wavelength-dependency). Because the angular in- 
formation of the 2007 observations is limited, we were not able 
to apply the same 2D model fitting strategy as above to those 
data (see Sect. [4]). Observed and modelled visibilities and differ- 
ential phases are displayed in the right panel of Fig. [2] We note 
a very good agreement between model and observations for the 
visibilities and for the differential phases. The model flux ratios 
of the Gaussian envelope and of the blob over the central star are 
plotted in Fig. [4] For comparison we show the contributions of 
relevant molecular and dust species that charac terise the spec- 
trum of a C-rich AGB star at these wavelengths (Now otnv et all 
201 1). The flux ratio of the blob over the central star follows the 
shape of the C2H2 contribution, with a slight shift of the mini- 
mum towards the shorter wavelengths (where the HCN contri- 
bution has a minimum). A weak drop can be recognised around 
1 1.3 fim where the SiC feature is located. The spectral feature in 
the differential phase (Fig. |2]i and the flux ratio (Fig.|4| suggest 
that the asymmetry originates in the C2H2+HCN layer (8-9 //m 
region), and possibly also in the SiC-forming region. The fea- 
tures of Fbiob/^star could be caused by the star or by the nature 
of the blob. This will be discussed in Sect. [4] 
The best-fitting model (reduced y 1 - 0.29, 4 free parameters) 
yields a FWHM for the Gaussian envelope of 70- 1 30 mas. These 
values are cons istent with single-di sh telescope mid-infrared ob- 
servations by Lag adec et alJ d201 ll) that showed an unresolved 
source with an angular resolution of ~ 0.3 arcsec. The model 
also yields a separation of the point source of 42 mas at an angle 
of -36°. We derived a Zf-band diameter of 8 mas (unpublished 
AMBER observations, PI Ohnaka), i.e. the asymmetric struc- 
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Fig. 2. The shaded-region show the calibrated visibilities (upper panels) and differential phases (lower panels) versus wavelength for the selected 
observations of 2007 (left) and 2009 (right). The overplotted lines show the best-fitting model (details in Sect.[3]l. 



ture is located at about 5 stellar radi i (~30 AU assuming 680 
pc as distance; IWhitelock et a l.. 2006) inside the envelope of the 
star. The position of the Delta function has to be treated with 
caution, because it is derived by fitting data with limited angu- 
lar coverage. The fit presented here provides one example of a 
quantitative solution, but more complex morphologies cannot be 
excluded. 



4. Discussion 

As already mentioned, a differential phase signature, similar to 
that observed for R For, was de tected for two J-typ e carbon stars. 
IRAS 1800-3213, presented bv lDeroo etaP d2007l) . shows a sig- 
nature in the differential phase that is in the same wavel ength 
region (i.e. 8-9 /mi) as the one observed for R For. lDeroo et alj 
(120071) mentioned that long baselines probe the part in the at- 
mosphere where dust emission is strongly (spatially) resolved. 
The authors speculated that the differential phase signature is 
caused by an offset of the photocentre between the strongly re- 
solved dust region and the unresolved stellar photosphere. This 
might be explained by a circumbinary disc. The second object, 
BM Gem, presented by Ohn aka et al.l (120081) . shows a non-zero 
differential phase signature between 9-11 /mi. The authors in- 
terpreted this as due to the presence of a circum-companion disc. 

The detected asymmetry for R For excludes the symmetric 
scenarios 1. and 2. hypothesised bv lFeast et al.1 (fl984). The plau- 
sible explanations are (i) a feature intrinsic to the star, such as a 
dust blob that moves through the circumstellar environment of 
R For; or (ii) an unresolved companion that is embedded in the 
circumstellar envelope. 

The dust blob might have been ejected in a preferential di- 
rection or randomly, as stated in scenarios 3 and 4. The ex- 
pansion velocity of the outer envelope of R For is 16.9 km s _1 



dGroenewegen et aTlll999l) . Using the determined separation of 
42 mas, and assuming that the blob moved with a constant ve- 
locity of 5kms _I over the last 30 years, (a reasonable value 
for th e inner envelope, pred icted by dynamic model computa- 
tions; iNowotnv et al.L HblO), it is worthwhile to mention that 
the event of 1983 would correspond to the moment at which 
such a structure was formed near the surfac e of the star. Recen t 
results from investigations of R CrB stars (Jeff ers et al.L 1201 2|) 
indicate the presence of dust clo uds with different grain-sizes. 
Nevertheless, Brig ht et al.l d201 lb did not detect any differential 
phase signature for these objects, arguing that the signatures of 
dust clumps are not strong enough to induce significant spec- 
tral signatures like the one we observed. This might not be the 
case for AGB stars because R CrB stars are much hotter and 
H-deficient (i.e. no C?H?), therefore ca rbon nucleation proceeds 
differently dGoeres & Sedlmavrl fl992b . IRC+10216 could be a 
very interesting test-case to check if the clumps can produce a 
phase signature in carbon stars. Indeed, several dust clumps were 
detected in th e atmo s phere of this star dHaniff & Buschen, 119981 ; 
IWeigelt et all I1998L bool iLeao et all 120061) at spatial scales 
comparable to the one we obtained with the GEM-FIND mod- 
elling (~ 30 AU). 

The binary option is quite unrealistic for a stellar compan- 
ion for the following reasons. Archive GALEX observations do 
not show any evidence of UV excess, which automatically ex- 
cludes the presence of a warm white dwarf. By translating the 
flux ratio (see Fig. |4]i into magnitudes, this would imply that 
the companion has a brightness of ~ 2.2 mag at 12 /mi. A faint 
giant star is very unlikely because such an object would be re- 
solved. A main sequence star also turns out to be too massive 
and thus would disrupt (because of mass accretion) the spher- 
ical symmetry detected at low spatial frequencies for the en- 
velope. We cannot exclude that the companion is a substellar 
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object orbiting the envelope. There are suggestions in litera- 
ture that such a s ituation might be common among AGB stars 
dSoszvnsku 120071). In this case, only low-mass accretion is ex- 
pected dLecave lier Des Etangs[ 120071) . leaving a signature only 
at the high spatial frequencies probed by long baselines. 

The differential phase signature was detected with the same 
baselines at the same position angles after two years. Assuming 
that the flux ratio of the point source over the central star is con- 
stant between 2007 and 2009 and that the position of the com- 
panion as well as the dust blob would change from 2007 to 2009, 
the following scenarios are possible: (i) for 1 Mq central star 
and a Jupiter mass companion (10~ 3 M G ) the orbital period of 
the companion would be 185 yr and would result in a change of 
the position angle of +4 °; (ii) assuming a dust blob that is mov- 
ing outwards at a uniform velocity of 5 km s _1 would result in a 
change of the separation of 2.4 mas. In both cases the asymmetry 
would be still observable within the two years with the chosen 
MIDI configuration. Therefore, the available observations do not 
allow us to reject any of these two asymmetric scenarios. 

This letter presents one interpretation for the detected asym- 
metry in the dust envelope of R For. More complex morpholo- 
gies cannot be excluded a priori because of the limited cover- 
age of the available observations. Interferometric imaging with 
the se cond-generation instrument VLTI/MATISSE dLopez et al.L 
2006) and a detailed polarimetric investigation are needed to 
clarify the nature of the detected structure. 
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Table 2. List of the calibrators. 



Calibrator 


Fn 


#UD 






[Jy] 


[mas] 




HD 224935 


86.90 


7.24 4 


0.03° 


HD 25025 


109.6 


8.76 + 


0.624 * 


HD 20720 


112.8 


10.03 


±0.1 


HD 12929 


77.80 


7.43 i 


0.52" 


HD 48915 


143.1 


5.85± 


0.15 a 



Table 3. Parametric description of the geometrical model 
CircGauss+UD+Dirac. 



Notes. The 12 /im flux listed in col. 2 is the IRAS flux. (n) ESO/MIDI 
database. JMMC database. M lDvck et al.Ul996h . 



Parameters 


Results 


slar , mas 


15-22 


FWHM, mas 


70-130 


^env /^star 


0.75-1.08 


-fblob/^star 


0.02-0.07 


p, mas 


42 





-36° 



100 



50 - 



-50 



-100 



1 1 1 1 1 1 1 1 1 




O' : 


. . . . i . . . . 





-100 -50 50 

u foresee"'] 



100 



Notes. A range of values is listed for the wavelength dependent pa- 
rameters. 9 ce „. .. diameter of central star; FWHM. ..full width at half 
maximum of the Gaussian envelope; Fj/Fj. . .flux ratios for Gaussian 
envelope/central star and blob/central star; p, (j>. .. separation and posi- 
ti on angle of the comp onent from the primary star. Detailed description 



tion angle ot the comp' 
in lKlotzetai] (l2012b). 



Fig. 1. uv coverage dispersed in wavelength of the observations ac- 
quired during the programs 080.D-0231 (grey-scale), and 084.D-0361 
(red-scale). The colour code ranges from 8-12.5 yum, with light colours 
indicating shorter wavelengths. The sample is restricted to the observa- 
tions used for this work (marked with (b) and (d) in TableQ}. 




Fig. 3. Normalised intensity distribution at 10 /im in the sky plane of 
the geometric model that best fits the R For MIDI data of 2009. The 
position of the asymmetric structure is indicated with a circle. The axis 
orientation is chosen according to the standard with north up and east 
left. 
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Table 1. Journal of the MIDI observations for the C-star R For. 



# 


Date & Time 


Phase* 


VLTI Configuration 


S P 


PA 


Seeing 




[UT] 






[m] 


[deg] 


[arcsec] 


01° 


2007-10-04 03:17:13 


1.07 (a) 


G0-H0 


26 


43 


0.82 


02" 


2007-10-05 04:58:33 


1.08 (a) 


G0-H0 


31 


61 


0.88 


03" 


2007-10-29 03:17:27 


1.14 (a) 


A0-K0 


122 


60 


1.13 


04 


2007-12-02 01:47:31 


1.23 (b) 


D0-H0 


63 


66 


2.05 


05 


2007-12-02 02:51:00 


1.23 (b) 


D0-H0 


63 


73 




06" 


2007-12-02 03:00:50 


1.23 (b) 


D0-H0 


63 


74 


_ 


07 


2007-12-06 01:10:17 


1.24(b) 


A0-K0 


125 


63 


1.35 


08 


2007-12-06 02:36:01 


1.24(b) 


A0-K0 


127 


73 


0.88 


09 


2007-12-10 00:58:46 


1.25 (b) 


DO-HO 


63 


64 


1.17 


10 


2007-12-10 01:50:54 


1.25 (b) 


DO-HO 


64 


70 


1.28 


11 


2007-12-12 03:11:41 


1.25 (b) 


DO-HO 


61 


79 


1.44 


12 


2007-12-29 01:54:08 


1.30(b) 


DO-HO 


61 


78 


1.32 


13 


2008-02-27 01:07:55 


1.45 (c) 


DO-HO 


67 


99 


0.99 


14" 


2009-09-06 07:32:39 


2.89 (d) 


DO-HO 


63 


67 


1.08 


15 


2009-09-06 09:06:52 


2.89 (d) 


DO-HO 


62 


77 


0.95 


16 


2009-09-07 05:43:04 


2.89 (d) 


GO-HO 


28 


52 


0.76 


17 


2009-09-07 06:18:38 


2.89 (d) 


GO-HO 


30 


57 


0.48 


18 


2009-09-07 09:10:11 


2.89 (d) 


HO-KO 


31 


78 


1.30 


19 


2009-09-07 09:44:55 


2.89 (d) 


HO-KO 


30 


81 


1.45 


20 


2009-09-21 07:53:15 


2.93 (d) 


A0-G1 


126 


76 


1.08 


21 


2009-09-26 09:00:58 


2.94 (d) 


D0-G1 


68 


152 


0.75 


22 


2009-09-27 06:01:02 


2.94 (d) 


D0-G1 


71 


130 


0.68 


23 


2009-09-27 08:43:21 


2.94 (d) 


DO-HO 


57 


82 


0.65 


24 


2009-10-01 05:21:32 


2.95 (d) 


EO-GO 


15 


63 


1.15 


25 


2009-10-01 05:59:46 


2.95 (d) 


EO-GO 


16 


67 


1.09 


26 


2009-11-18 03:31:27 


3.08 (a') 


A0-G1 


128 


72 


0.90 


27 


2009-11-18 03:50:10 


3.08 (a') 


A0-G1 


127 


74 


0.65 


28 


2010-01-11 03:13:02 


3.22 (b') 


DO-HO 


38 


96 


0.65 


29 


2010-01-11 03:48:17 


3.22 (b') 


DO-HO 


45 


92 


0.76 



Notes. w Estimated from AAVSO light curve (http://www.aavso.org/aavso-research-portal). (a) Data exclude d from astrophysic al 
interpretation, because they were not fulfilling one or more of the five criteria for the quality assessment listed in Sect. 2.2 of Kl otz et all i2012ah . 



